Dynamics of long-range ordering in an exciton-polariton condensate 
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We report on time resolved measurements of the first order spatial coherence in an exciton- 
polariton Bose-Einstein condensate. Long range spatial coherence is found to set in right at the onset 
of stimulated scattering, on a picosecond time scale. The coherence reaches its maximum value after 
the population and decays slower, staying up to a few hundred of picoseconds. This behavior can 
be qualitatively reproduced, using a stochastic classical field model describing interaction between 
the polariton condensate and the exciton reservoir within a disordered potential. 
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Long range spatial coherence lays at the heart of the re- 
markable quantum behaviour of matter manifested by su- 
perconductors and superfiuids at the macroscopic scale. 
An intriguing question is how quickly the coherent field 
can be established when a normal incoherent gas is driven 
through the phase transition. Experimentally, such dy- 
namics are not easily accessible in traditional supercon- 
ductors and superfiuids. Only recently the question of 
the formation of long range coherence was addressed in 
experiments with ultracold atomic gases [HH], while the 
importance of this study was pointed out even before the 
realization of Bose Einstein condensation (BEG) |3] |J . 

The demonstration of BEG in a solid state system was 
achieved a few years ago with exciton-polaritons [5] . An 
exciton-polariton is a quasi-particle that appears in the 
regime of strong coupling between a microcavity photon 
and a quantum well exciton. It inherits desirable prop- 
erties for condensation from both components. Their ex- 
citonic content allows for the collisional relaxation, while 
the photonic component provides a light mass. 

Moreover, the polariton statistical properties are di- 
rectly accessible through the photoluminescence due to 
a linear relation between external and internal cavity 
fields [B]. The cavity radiation field is a part of the po- 
lariton wave function. Therefore, the spatial phase co- 
herence, which is the order parameter of a BEG tran- 
sition [7], can be directly probed. In our experimen- 
tal work we experimentally measure the first order cor- 
relation function in a direct way. Time resolved de- 
tection setup in combination with pulsed excitation al- 
lows us to observe how the spatial phase correlations 
are built in time. We therefore probe the time evolu- 
tion of the first order spatial coherence function [H [9]: 
g^^\r,r',T = 0,t) =< ip* {v,t)il}{v' ,t) >, at zero delay 
time r. 

Our sample is the same one as in our previous works [5 , 
[TUJ [TT]. It is a GdTe/GdMgTe microcavity grown by 
molecular beam epitaxy that contains 16 quantum wells 
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with a vacuum field Rabi splitting of 26meV. The sam- 
ple was mounted in a liquid-He flow cryostat and cooled 
down to approx. lOK. The sample was excited through 
a 50x microscope objective of numerical aperture NA = 
0.5 providing a diffraction-limited spatial resolution. The 
excitation spot was tuned to ~ 15/im diameter in a quasi 
top-hat intensity profile. We excited the sample non- 
resonantly with a pulsed Ti: Sapphire laser (250fs pulse 
duration, / = 80MHz repetition rate) tuned to the res- 
onance with energy of the first minimum in the Bragg 
mirrors at 695nm. In this way we create an electron - 
hole population above the semiconductor band gap. The 
photo-created electron - hole pairs first relax into excitons 
through scattering with optical phonons. These excitons 
form a thermalized reservoir with a long lifetime of hun- 
dreds of ps [T21[T3J[I1]. When the exciton density is suffi- 
ciently high, the exciton-exciton scattering rate into the 
LP branch overcomes the radiative losses and the pro- 
cess of stimulated scattering starts [13]. The polariton 
population lasts much longer than the polariton life time 
(few picoseconds). This is possible because the excitonic 
reservoir keeps replenishing the polariton population. 

The coherence of the lower polariton population is 
probed with a Michelson interferometer in the mirror- 
mirror configuration. We superimpose the polariton pho- 
toluminescence with its copy translated by a distance 
of Ax, by tilting one of the mirrors (Fig. [T|). In this 
way, we can investigate the spatial correlations between 
two points from spatially distant parts of the condensate. 
Below threshold, the system exhibits short-range order, 
with a correlation length of ^ 0.8/im. Ghoosing the dis- 
tance Ax — 8.5/im, we probe the coherence properties of 
points separated by more than ten times the thermal de 
Broglie wavelength. To trace in time the appearance of 
the condensed phase we have selected a slice in the in- 
terference pattern, as it is shown in Fig. [l] and mapped 
out its temporal evolution with a streak camera. To re- 
late our temporally resolved interferogram with the first 
order correlation function we scanned the phase in the 
Michelson interferometer over ~ Stt (see Ref. |5] for more 
details) . 
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FIG. 1; Time integrated real space image of the exciton- 
polariton condensate emission at the output of the Michel- 
son interferometer. The signal reflected by the two mirrors is 
marked by black circles. The doted vertical lines are guides to 
the eye coinciding with the minima of the interference fringes. 
The part of the photoluminescence in the middle of the over- 
lapped region, at the position of the white solid line, was sent 
to a streak camera and resolved temporally. 



Fig. |2] illustrates the time-resolved images of the 
buildup of coherence over the selected line in the con- 
densate. The appearance of the long range coherence is a 
manifestation of the phase transition. If the power in the 
excitation pulse is not high enough to provide sufficient 
number of polaritons to cross the condensation thresh- 
old, the coherent state is not formed (Fig. [2^). Above 
threshold, the coherence is formed after a certain delay 
from the excitation pulse and can persist for more than 
lOOps. The higher the excitation power, the shorter the 
delay of the coherence formation. The full dynamics of 
the coherence (buildup and decay) is strongly accelerated 
as we increase the excitation power. This acceleration is 
due to an increasing exciton-exciton collision rate and is 
an indirect proof of the fact that exciton-exciton colli- 
sions play an important role in the exciton relaxation. 

The inhomogeneous distribution of photoluminescence 
intensity within the excitation spot (Fig. [T|) and the co- 
herence in the spatial direction results from the disorder 
in the sample that modulates the polariton density [TU]. 
The time evolution of coherence shows also an irregular 
behaviour. The higher the excitation power the more 
the coherence fluctuates in time. This effect is due to re- 
distribution of polaritons over the excitation spot as the 
effective disorder seen by the particles is changing with 
the population [10]. To get a general understanding of 
the coherence behaviour, less influenced by the disorder, 
we have averaged the observed coherence over the differ- 
ent correlated points, that always have the same spatial 
separation of Ax = 8.5/im between them (see Fig. |3]). 

The time evolution of the polariton density is traced 



(Fig.jsjred line) simultaneously with the coherence (Fig.js] 
black dots), by temporally resolving the emission re- 
flected separately by the two mirrors of the Michelson 
interferometer. This gives a more detailed insight in the 
time at which the phase coherence starts to appear with 
respect to the particle population. The dashed vertical 
lines in Fig. [3] separate the low polariton density region, 
from the region of high density. At this border the den- 
sity evolution shows a cusp, indicating the onset of stim- 
ulated scattering. Exactly at the same time (below our 
temporal resolution limit of 3.5ps), the coherence starts 
to build up. This directly shows that the communication 
between points separated by 8.5/xm in the condensate so 
as to define a common phase, is very fast, in a picosecond 
time scale. It gives a phase communication velocity faster 
than 2.5 • lO^m/s. The total time span of the condensate 
formation can be evaluated by the time needed for an 
increase of the coherence from 10% to 90% of the maxi- 
mal value. We thus obtain 12ps at the threshold Pq, and 
6.5ps for an excitation power of 1.4-Po. This means that 
the total coherence is established here with a velocity 
0.7- lO^m/s at threshold and 1.3 • lO^m/s for 1.4-Po. The 
higher the excitation power the faster the coherent phase 
is built. These velocities are faster or comparable with 
the velocity of sound in the polariton BEC, c = ^/Jl/rn. 
We estimate c from the blue shift /i induced by the con- 
densate polaritons and the lower polariton mass to, to be 
(0.7 ± 0.2) • lO^m/s. This means that the appearance of 
coherent phase is faster than any process mediated by the 
interactions. Our result differs in this case from the cold 
atomic gases [T], where it was found that the coherent 
region grows with a velocity of approximately a factor 
of five slower than the speed of sound in atomic BEC. 
Atomic BEC are closed systems, therefore the only way 
for correlations to propagate is using interactions within 
the atomic cloud. Polaritons form a driven dissipative 
system. Moreover, the condensed fraction of polaritons 
is in constant interactions with the excitonic reservoir 
and phonons. Therefore the correlations can be indeed 
established faster than the speed of sound. Once it is 
established, the phase coherence slowly decays, and re- 
markably enough, it is maintained in the system for a 
long time, of about one hundred ps, decaying slower than 
the population for high excitation powers. 

A comparison of the data measured at the different 
pump powers can be made in Fig. [4] where the coherence 
is plotted versus the density for the three different pump 
intensities. This figure shows the importance of the speed 
at which the phase transition is crossed. At low pump in- 
tensity (green curve) , where the density builds up slowly, 
the coherence is almost a single valued function of the 
condensate density. Physically, this can be interpreted 
as an adiabatic evolution of the coherence with the den- 
sity. At higher powers (blue curve) on the other hand, a 
hysteresis loop opens up. At initial times, the coherence 
is below the steady state value. During the fast density 
ramp, the phase coherence does not have time to spread 
through the system. 



3 



y spatial direction Ipm) 
6 12-12 -6 




FIG. 2: Dynamics of first or- 
der correlation function, g'^', be- 
tween two lines of the conden- 
sate separated by 8.5/im for exci- 
tation powers (a) 0.5-Po, (b) Po, 
(c) 1.4-Po and (d) 2-Po, where 
Po is the threshold power. The 
colour scale denotes the g^^^ am- 
plitude, white is used at the po- 
sitions where the g'^' can not be 
defined or at the position of the 
laser pulse at time t =Ops. 



The same qualitative behaviour is obtained in theo- 
retical simulations W? . Our model describes the lower 
polariton dynamics by a stochastic classical field equation 
coupled to a rate equation for the exciton reservoir. The 
latter is split in two parts: the active excitons that can 
match energy and momentum conservation for a scat- 
tering process into the lower polariton branch (density 
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FIG. 3: Transient of population and coherence, solid red line 
and black dots, respectively, for the average over along 
the white solid line in Figjll for excitation powers (a) 0.5-Po, 
(b) Po, (c) 1.4-Po and (d) 2-Po, where Po is the threshold 
power. The grey line shows the estimated standard deviation 
of the fit coefficients obtained in the g'^' calculation, (e) and 
(f) illustrate the theoretical simulation results at and above 
threshold, respectively. Depicted here is the total polariton 
population at —{Ax/2) and {Ax/2) positions. The onset of 
the stimulated scattering is marked by the vertical dashed 
lines. 



nR2) and the excitons that cannot directly scatter into 
the lower polariton branch (density n^i). The coupled 
equations of motion of the reservoirs describe the non- 
resonant excitation with intensity P, the equilibration of 
the densities between the two reservoirs at a rate K12 and 
the scattering into the lower polariton region 



dt 
d 
dt 



nm = P - Ki2{ nR2)-lRinRi 

I nRi , dnR2 

nR2 = Ki2{ nR2) - jR2nR2 - a 

ri2 dt 



(1) 
(2) 

R2^LP 



where ri2 is the equilibrium ratio between the two 
excitonic reservoirs and the last term on the sec- 
ond line accounts for the scattering into the lower 
polariton region. The coupling terms between the 
reservoir and lower polariton branch read explicitely 

'7^in,out'/'(x) = e""'''V"-fl2(x)nfl2(x')-Rin,out(eq)V'x' 

The second reservoir density is normalized so that thresh- 
old is reached when it is unity. The parameter a quanti- 
fies the depletion of the exciton reservoir due to its scat- 
tering into the lower polariton field -0, defined on a spa- 
tial grid with unit cell area AV^, that is described by the 
stochastic classical field equation 



id^ = dt 



7) 



2m 



Kxt5lV'(x)|' 



+dW, 



where Voxt is the disorder potential that acts on the po- 
laritons. The noise term dW is a complex Gaussian 
stochastic variable with the correlation functions propor- 
tional to the sum of loss and gain. The advantage of this 
method for the modeling of the experiments is that spa- 
tially inhomogeneous systems can be described without 
any additional cost. We are thus able to model a dis- 
ordered microcavity excited with a spatially finite laser 
pulse. 

Figs. [Sjj, f show the results for two different excitation 
intensities |16j . They show qualitatively the same trend 
as the experimental data. Close to threshold (Fig. [3^), 
the coherence follows the population, while it stays be- 
hind further above threshold (Fig. |3|^), and is kept for 
longer times. In the coherence versus population graph 
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FIG. 4: (a) Experimental and (b) theoretical results of build 
up and decay of coherence as the population of polaritons is 
formed at threshold (green) , above threshold (blue) and twice 
above threshold (red). The timeline is marked by arrows. 
For high excitation powers we evidence the appearance of a 
hysteresis between the rise and decay. 



(Fig. |4p), the hysteresis loops appear (blue curve) and 
as in the experimental data its upper part overlaps with 
the curve at lower pump intensity (green curve). 

The disorder potential limits the coherence of the po- 
lariton condensate. Simulations without disorder predict 
a maximal coherence of about 0.45. Less expected is that 
the finite size excitation spot also plays a role in the co- 
herence formation time. With respect to an infinite sys- 
tem (simulated by using periodic boundary conditions), 
the coherence builds up faster in the finite system. The 
physical reason is that spontaneously formed vortices in 
the condensate phase [T71 |T51 HH] can escape away from 
the system at its boundaries, instead of only annihilat- 
ing each other. Finally, simulations have shown that the 
out-scattering term in the model (see equation ^) plays 
an important twofold role. On one hand, it decreases the 
time needed to build up the coherence, on the other hand, 
it decreases the maximal value of the coherence that is 
achieved [15] . The latter effect is due to the enhancement 
of fluctuations, whereas the initial acceleration of the co- 
herence formation is due to the out-scattering contribu- 
tion to the deterministic term that quickly eliminates the 
high momentum components of the polariton gas. Con- 
ceptually, this is important, because it shows that the 
coherence of the low energy polaritons is strongly influ- 
enced by their interaction with the excitonic reservoir. 
The polariton subsystem behaves in this way differently 
from ultracold atomic gases, which are isolated from their 
environment. 



Let us now comment on the experimental data at the 
highest pump intensity (Fig. [sji and the red curve in 
Fig. [4^). In this ultrafast regime, the coherence builds up 
surprisingly fast. Not only does the hysteris loop (Fig.|4|i 
red curve) turn in the other direction with respect to the 
blue curve, but the coherence at short times is even above 
the value that is reached in the adiabatic regime. Within 
our theoretical model, we have not been able to find a 
regime where this behaviour is realized. It is possible 
that it originates from a nontrivial dynamics in the exci- 
ton reservoir, but could also be related to the finite time 
measurements. Note that the error on the coherence at 
short times (gray line in Fig. [sji) is large at the initial 
time where the coherence is formed. 

In conclusion, our work has provided substantially new 
experimental observations on the formation dynamics of 
the long range order in the exciton-polariton BEC. The 
all-optical access of the polariton condensate properties 
is a clear advantage to study the dynamics of its coher- 
ence formation. We have demonstrated that the onset 
of phase coherence between points separated by 8.5/im 
in the condensate appears in a picosecond time scale, 
much faster than expected: simultaneously with the on- 
set of stimulated scattering. When the dynamics in the 
system is slow, at threshold, the coherence adiabatically 
follows the particle population. When the system is 
rapidly driven through the phase transition, the maxi- 
mum long-range order is established in the system slower 
than the buildup of the population. This demonstrates 
that the system needs some time, a few picoseconds in 
our case, to buildup the highest long range order. Once 
it is formed, the coherence can be maintained in the po- 
lariton system for more than a hundred picoseconds. We 
have found that the velocity of the establishment of the 
phase in a polariton condensate is as fast or faster than 
the sound velocity, i.e. much faster than in cold atom 
systems. These experimental results are supported by a 
theoretical model which brings new elements to the un- 
derstanding of the mechanisms that regulate the phase 
mediation in non-equilibrium condensates. The observa- 
tion and the understanding of these ultrafast dynamics 
are also of crucial importance for the design of future 
polariton laser devices. 
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